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A  COMPUTER  PROGRAM  FOR  INTERPOLATION 
AND  PLOTTING  OF  ISOHYETAL  POINTS^ 

M.  H.  Disking  and  P.  R.  Davis^ 


INTRODUCTION 


An  important  part  of  the  study  and  interpretation 
of  rainfall  records  is  the  preparation  of  isohyetal  maps. 
Such  maps,  showing  hnes  of  equal  depth  of  rainfall,  may 
be  prepared  for  the  total  storm  rainfall,  for  rainfall 
during  selected  time  intervals  within  the  storm,  or  for 
rainfall  accumulated  from  a  number  of  storms  over  a 
specified  period  of  time  such  as  a  month,  a  season,  or  a 
year. 

Manual  preparation  of  isohyetal  maps  is  tedious, 
time  consuming,  and  subjective.  The  subjectivity  is 
largely  due  to  the  interpretation  of  the  original  data  by 
the  person  preparing  the  map.  This  is  not  entirely  a 
disadvantage,  since  the  number  of  rain-gaging  stations  is 
usually  less  than  desirable  for  preparation  of  good 
isohyetal  maps,  and  some  human  judgment  is  needed  for 
the  production  of  acceptable  or  reasonable  maps,  based 
on  the  insufficient  data. 

Some  computer  programs  for  the  preparation  and 
plotting  of  isohyetal  maps  are  commercially  available. 
The  programs  are,  however,  expensive  both  in  terms  of 
the  initial  investment  needed  to  develop  or  to  procure 
the  program  and  in  terms  of  the  computer  time  needed 
for  the  preparation  of  an  isohyetal  map.  One  reason  for 
the  complexity  and  slowness  of  these  programs  is  the 
need  to  provide  for  all  possible  combinations  of  deci- 
sions that  make  up  the  factor  referred  to  as  human 
judgment. 

The  compromise  between  the  manual  plotting  of 
isohyetal  maps  and  the  computer  preparation  of  such 
maps  is  a  computer  program  that  carries  out  the  tedious 
and  mechanical  stage  of  the  map  preparation,  leaving  to 
the  investigator  or  draftsman  the  stage  in  which  human 
judgment  is  to  an  advantage.  The  first  of  these  two 
stages  consists  of  the  interpolation  between  rain-gaging 
stations  to  find  points  where  the  rainfall  has  values  equal 
to  integral  multiples  of  the  contour  interval  selected  for 
the  map  and  the  plotting  of  these  isohyetal  points.  The 


second  stage  involves  the  plotting  of  the  isohyets  in  the 
form  of  smooth  curves  that  pass  through  or  near  the 
isohyetal  points  determined  in  the  first  stage. 

Assigning  the  interpolation  and  plotting  of  iso- 
hyetal points  to  a  computer  program  also  contributes 
appreciably  toward  a  reduction  in  the  subjectivity  of  the 
resulting  isohyetal  maps.  Experience  with  the  prepara- 
tion of  such  maps  indicates  that  many  of  the  differences 
between  maps  produced  by  different  individuals  can  be 
traced  to  the  way  they  prepare  and  plot  the  interpolated 
isohyetal  points.  Since  this  part  of  the  process  is  tedious 
and  uninspiring,  many  individuals  tend  to  estimate 
rather  than  compute  the  location  of  the  points  and  to 
reduce  the  number  of  such  points  to  a  bare  minimum. 
The  resulting  map  reflects,  of  course,  the  decisions  made 
at  this  stage.  A  computer  doing  this  part  of  the  work 
increases  the  accuracy  of  the  plotted  points  and  makes  it 
possible  to  plot  all  the  points  needed  for  the  preparation 
of  the  isohyetal  map.  It  is  even  possible  to  produce  more 
points  than  are  actually  needed  for  the  tracing  of  the 
isohyets.  The  extra  points  are  obtained  by  specifying  in 
the  program  a  smaller  contour  interval  and  additional 
lines  of  interpolation.  They  are  useful  in  making  the 
decisions  about  the  most  probable  configuration  of  the 
isohyets. 

A  program  for  the  interpolation  of  isohyetal 
points  between  rain  gages  and  for  the  plotting  of  these 
points  at  a  desired  scale  has  been  developed  jointly  at 
the  Southwest  Watershed  Research  Center,  Agricultural 
Research  Service,  and  the  Department  of  Agricultural 
Engineering,  University  of  Arizona.  The  program  is 
based  on  two  previous  programs  for  the  preparation  of 
ground-water-level  contour  maps  developed  by  the  De- 
partment of  Agricultural  Engineering,  University  of 
Arizona.  A  description  of  the  new  program,  some 
examples  of  its  use,  and  information  about  the  computer 
time  needed  for  its  operation  are  discussed. 
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THE  INTERPOLATION  PROCEDURE 


The  interpolation  procedure  adopted  for  the  com- 
puter program  is  the  same  as  that  used  in  manual 
preparation  of  isohyetal  maps.  It  consists  of  a  simple 
linear  interpolation  along  straight  lines  connecting  adja- 
cent stations.  Thus,  if  the  depths  of  rainfall  at  any  two 
adjacent  stations  are  Hj  and  H2  ,  and  the  distances  from 
the  stations  to  a  point  on  the  line  between  them  are  Li 
and  L2  ,  the  interpolated  depth  of  rainfall  at  the  point  is 

H  =  Hi+,     1,  (H2-H,)  (1) 


or 


H  = 


Li  +  L2 


L2  Hi  +  Li  H2 
L,   +  L, 


(2) 


Denoting  the  distance  between  the  two  stations  by  L, 
the  equations  may  be  solved  for  the  distance  along  the 
line  from  Station  No.  1  to  a  point  where  the  inter- 
polated depth  is  a  given  value  H. 


L,  = 


H  -  Hi 
H2  -Hi 


(3) 


If  the  stations  are  specified  in  terms  of  their 
coordinates  on  a  rectangular  X-Y  axes  system  (fig.  1), 
the  equations  may  be  rewritten  in  terms  of  these 
coordinates.  The  equations  for  the  coordinates  (Xp,  Yp) 
of  the  point  where  the  depth  of  rainfall  is  H,  correspond- 
ing to  equation  3,  are 


and 


_    H-Hi 

Xp  -  Xi (Yi  -  Xi ) 

H,  -H, 


H-Hi 
Yi  =  „      „      (Y2-Y1) 


(4) 


H2  -Hi 


(5) 


The  equations  are,  of  course,  not  restricted  to 
points  between  the  two  rain-gaging  stations.  They  may 
be  used  also  for  extrapolating  to  points  situated  on  the 
straight  line  joining  the  stations  but  outside  the  intercept 
defined  by  the  two  stations. 

The  first  step  in  the  application  of  the  above 
interpolation  procedure  to  the  preparation  of  an  iso- 
hyetal map  is  the  construction  of  a  network  of  lines 
connecting  adjacent  stations.  This  network,  which  may 
be  called  the  interpolation  network,  is  composed  of  an 
assembly  of  lines  along  which  the  interpolation  proce- 
dure is  applied.  Ideally,  the  interpolation  network  would 
be  composed  of  an  assembly  of  equilateral  triangles  (fig. 
2a),  which  could  be  obtained  if  the  rain-gaging  stations 
were   situated  on  a   triangular  grid.   In   practice,  this 


(0) 


Figure  1. -Interpolation  for  a  line  joining  two  stations. 


(b) 

Figure  2. -Various  shapes  of  interpolation  networks:  (a)  equi- 
laterial  triangles,  (b)  four-sided  polygons,  and  (c)  five-sided 
polygons. 


configuration  is  very  rare,  and  the  network  will  usually 
be  composed  of  four-  or  five-sided  polygons,  for  which 
interpolation  may  be  carried  out  along  two  or  more 
intersecting  diagonals  (figs.  2b  and  2c). 

Interpolated  values  of  the  depths  of  rainfall  at  the 
point  of  intersection  of  any  two  diagonals  will  usually 
not  be  equal,  depending  on  the  values  of  rainfall  at  the 
stations  joined  by  these  diagonals.  In  the  manual 
preparation  of  isohyetal  maps,  one  line  (or  at  least  a  part 
of  the  Une)  of  each  pair  of  intersecting  diagonals  is 
frequently  ignored.  The  decision  of  which  diagonal  to 


ignore  is  often  arbitrary  and  may  lead  to  errors  in  the 
final  map.  With  a  computer  interpolation  program,  it  is 
possible,  however,  to  include  all  diagonals  and  to  obtain 
interpolated  values  along  all  such  lines.  The  decision  to 
reject  any  interpolated  point  or  a  group  of  such  points 
can  be  made  at  a  later  stage  when  all  the  points  are 
available  for  inspection.  It  should  be  noted  that  in 
different  rainfall  events  the  rejected  points  will  be 
located  in  different  parts  of  the  interpolation  network, 
depending  on  the  pattern  of  rainfall  recorded  at  the 
rain-gaging  stations. 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM 


The  computer  program  developed  is  made  up  of 
three  major  divisions:  (1)  The  first  examines  the  data 
and  computes  the  coordinates  of  interpolated  points 
along  hnes  of  the  interpolation  network  as  specified  in 
the  input  data;  (2)  the  second  rearranges  the  interpo- 
lated points  and  the  rain-gaging  station  in  the  order  of 
the  magnitude  of  their  X  or  Y  coordinates  as  a 
preparation  for  the  plotting;  (3)  the  third  plots  the 
rain-gaging  stations  and  the  interpolated  points  in  their 
locations  relative  to  a  rectangular  grid  at  the  scale  chosen 
for  the  map.  More  details  about  the  three  divisions  are 
given  in  the  following  paragraphs.  A  flow  chart  of  the 
computer  program  is  given  in  figure  3. 
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Figure  3.-Flow  chart  for  computer  program. 


In  the  first  division,  the  input  data  are  examined 
to  determine  if  the  coordinates  of  any  station  are  given 
more  than  once  or  if  there  is  more  than  one  value  of 
rainfall  per  storm  for  any  of  the  gages.  If  either  of  these 
cases  occurs,  the  program  will  print  a  note  that  there  are 
duphcate  data  and  will  stop.  During  the  preliminary 
examination  of  data,  the  program  also  detects  the 
presence  of:  (1)  gages  for  which  coordinates  are 
available  with  no  rainfall  data,  or  (2)  gages  for  which  the 
rainfall  is  given  but  no  coordinates  are  available.  In 
either  case,  a  note  will  be  printed  giving  the  number  of 
the  gage  affected,  but  the  program  will  continue  with 
the  computations  for  the  gages  for  which  the  data  are 
complete.  Stations  with  zero  rainfall  are  treated  the 
same  as  any  other  station.  During  execution  of  the  first 
division,  notes  will  also  be  printed  giving  the  number  of 
any  gage  for  which  interpolation  cannot  be  completed 
because  of  either  of  the  above  reasons. 

After  examination  of  the  input  data,  the  program 
computes  the  coordinates  of  the  interpolated  isohyetal 
points.  The  depths  of  rainfall  for  the  isohyetal  points  are 
taken  at  integer  multiples  of  the  contour  interval 
specified  for  the  map,  and  the  coordinates  are  computed 
by  equations  4  and  5.  Interpolation  is  done  along  the 
Unes  specified,  proceeding  from  the  station  with  lower 
rainfall  to  the  station  with  the  higher  rainfall.  Isohyetal 
points  falling  outside  the  segment  defined  by  the  two 
stations  are  also  included,  to  avoid  loss  of  information,  if 
their  distance  from  the  stations  is  less  than  3  percent  of 
the  distance  between  the  stations  or  two  plotting 
increments,  whichever  is  longer.  No  interpolation  will  be 
carried  out  between  two  stations  for  which  the  absolute 
value  of  the  difference  in  depth  of  rainfall  is  less  than 
0.001  of  the  contour  interval.  The  interpolation  proce- 
dure is  applied  to  each  station  on  the  watershed  in  turn. 
Each  of  the  stations  can  have  up  to  four  adjacent 
stations  to  which  the  interpolation  will  be  carried  out. 

In  the  second  division,  the  program  considers  all 
the  isohyetal  points  generated  in  the  first  part,  as  well  as 


all  gaging  stations  used  in  the  interpolation  process,  and 
prepares  them  for  the  third  part  of  tlie  program.  All 
points  in  these  two  groups  are  rearranged  in  a  decreasing 
order  of  the  value  of  their  Y  coordinate.  If  a  rotation 
option  has  been  specified,  whereby  the  map  is  to  be 
rotated  by  90°  before  printing,  the  rearrangement  is 
done  with  respect  to  an  increasing  value  of  the  X 
coordinate. 

The  list  of  stations  used  in  the  interpolation 
procedure,  as  well  as  their  coordinates,  and  the  values  of 
the  depth  of  rainfall  for  the  storm  considered  are  printed 
during  the  second  division  of  the  program.  Another  list 
of  interpolated  points  rearranged  as  described  above,  can 
also  be  printed  at  this  stage  if  a  print  option  is  specified. 

In  the  third  division,  the  program  prepares  and 
plots  the  map  consisting  of  rain-gaging  stations  and 
interpolated  isohyetal  points  plotted  at  their  proper 
locations  according  to  the  scale  of  the  map  used.  An 
example  of  the  map  produced  is  given  as  figure  7  in  the 
section.  Applications  of  the  Program.  The  rain-gaging 
stations  are  represented  by  asterisks  (*),  and  the 
isohyetal  points  are  represented  by  letters  of  the 
alphabet  (A,  B,  C,  D,  etc.).  The  letter  printed  at  each 
point  represents  the  value  of  the  isohyet  passing  through 
that  point.  Thus,  the  letter  A  represents  points  for  which 
the  interpolated  depth  of  rainfall  is  equal  to  one 
isohyetal  interval,  the  letter  B  represents  points  for 
which  the  depth  of  rainfall  is  equal  to  two  isohyetal 
intervals,  etc.  Points  with  zero  rainfall  are  indicated  by 
the  letter  "Z."  The  same  letter  is  used  also  for  points 
where   the   depth   of  rainfall  is  equal  to  26  isohyetal 


intervals.  If  the  range  of  rainfall  values  to  be  represented 
is  more  than  26  isohyetal  intervals,  the  sequence  of 
letters  repeats  itself. 

The  map  is  plotted  at  the  scale  specified,  with  the 
X  axis  across  the  printed  page  of  output  and  values  of  X 
increasing  from  left  to  right.  The  Y  axis  is  down  the 
printed  page  and  is  continuous  from  one  page  to  the 
next.  Values  of  Y  increase  from  the  bottom  of  the  page 
upwards.  If  the  range  of  X  values  exceeds  the  width  of 
the  printed  page,  additional  strips  of  the  map  will  be 
produced  so  that  the  full  range  of  X  values  is  covered. 
An  option  is  available  in  the  program  to  rotate  the  axes 
clockwise  by  90°  so  that  the  Y  axis  runs  across  the  page 
and  the  X  axis  is  perpendicular  to  the  width  of  the 
printed  page. 

Identification  of  the  watershed,  the  storm  event, 
the  scale  of  the  map,  and  the  contour  interval  used  are 
printed  below  the  map.  The  axis  parallel  to  the  width  of 
the  paper  is  labeled  as  the  X  axis  or,  if  the  map  has  been 
rotated,  as  the  Y  axis.  Both  axes  have  the  scale  in  miles 
printed  along  the  borders  of  the  map.  The  axis  across  the 
page  has  values  at  1 .0-inch  intervals.  The  axis  perpendic- 
ular to  the  width  has  the  coordinate  given  for  each 
printed  line.  The  value  given  is  the  value  of  the 
coordinate  at  the  center  of  the  printed  line.  The  number 
of  miles  represented  by  1  inch  is  the  same  along  the  two 
axes,  so  that  the  map  produced  is  free  of  distortion. 

A  listing  of  the  program  is  given  in  the  Appendix. 
Copies  of  the  deck  of  cards  of  the  program  are  available 
from  the  Southwest  Watershed  Research  Center,  Agricul- 
tural Research  Service,  Tucson,  Ariz.  85705. 


CAPACITY  AND  ACCURACY  OF  PROGRAM 


The  computer  program  in  its  present  form  can 
prepare  a  map  of  isohyetal  points  for  networks  com- 
posed of  up  to  200  rain-gaging  stations.  Each  of  these 
stations  can  be  the  origin  of  up  to  four  lines  to  adjacent 
stations,  along  which  interpolation  is  carried  out.  Since 
each  of  the  lines  in  the  interpolation  network  needs  to 
be  defined  only  once,  it  follows  that  the  total  number  of 
interpolation  lines  that  originate  or  end  at  a  given  station 
can  be  as  high  as  eight.  Experience  with  the  program 
indicated  that  this  number  is  quite  adequate. 

The  number  of  interpolated  isohyetal  points  gener- 
ated by  the  program  can  be  up  to  3,000  minus  the 
number  of  rain-gaging  stations  for  which  coordinates  are 
given.  This  number  allows  an  average  number  of  at  least 
14  interpolated  points  for  each  rain-gaging  station.  This 
appears  to  be  more  than  adequate. 

Both  limits  given  above  are  controlled  by  "DI- 
MENSION" and  "DO"  statements  within  the  program. 


The  two  limits  may  be  changed  without  impairing  the 
performance  of  the  program.  The  program  length  with 
the  above  limits  is  about  62,200  (Octal)  including  input 
and  output  buffers. 

The  accuracy  of  the  plotting  of  the  locations  of 
interpolated  points  and  of  the  rain-gaging  stations 
depends  on  the  scale  chosen  for  the  map.  The  plotted 
points  are  represented  by  printed  letters  on  the  output 
printer  of  the  computer.  The  maximum  error  possible  in 
the  true  location  of  the  interpolated  points  in  a  direction 
parallel  to  the  width  of  the  paper  is  therefore  half  the 
spacing  of  the  letters.  The  maximum  error  of  plotting  in 
the  direction  perpendicular  to  the  width  is  equal  to  half 
the  spacing  between  the  centers  of  printed  lines. 

The  program  was  written  for  a  printer  producing  6 
lines  per  inch  down  the  page  and  10  characters  per  inch 
of  printed  line.  Thus,  the  area  represented  by  any  one 
letter  has  the  dimensions 


and 


DELX  =  SCALE/ 10.0 


DELY  =  SCALE/6.0 


(6) 


(7) 


across  and  down  the  page,  respectively.  In  these  equa- 
tions, "SCALE"  stands  for  the  number  of  miles 
represented  by  1  inch  on  the  map.  If  the  program  is  run 
on  a  computer  having  different  spacings,  the  program 
instructions  corresponding  to  equations  6  and  7  must  be 
changed  to  agree  with  the  spacings  of  lines  and  letters  in 
the  computer  used.  If  the  number  of  letters  is  other  than 
10  characters  per  inch,  changes  will  also  be  necessary  in 
the  part  of  the  program  that  arranges  for  the  printing  of 
the  upper  and  lower  borders  of  the  map  and  the  scales 
printed  with  them. 

The  computed  locations  of  two  or  more  isohyetal 
points  may  be  such  that  the  letters  representing  them 
should  be  printed  at  the  same  location.  To  avoid  loss  of 
information  resulting  from  the  printing  of  only  one  of 
the  letters,  provision  is  made  in  the  program  to  print  the 
second  character  one  plotting  increment  to  the  riglit  of 


the  position  where  it  should  have  been  printed.  If  the 
two  locations  are  occupied,  the  point  will  be  plotted  to 
the  left  of  the  original  location.  The  maximum  plotting 
error  for  a  letter  that  is  part  of  a  group  of  adjacent 
letters  may  thus  be  as  much  as  1 .5  plotting  increments 
from  its  true  position. 

The  accuracy  of  plotted  isohyetal  points  can 
obviously  be  increased  by  using  a  larger  scale  (smaller 
number  of  miles  per  inch).  The  accuracy  of  drawing  the 
isohyetal  lines  depends  on  the  number  of  interpolated 
points  available.  It  was  found  that  the  accuracy  of  the 
isohyets  can  be  improved  by  specifying  a  contour 
interval  equal  to  half  the  contour  interval  needed  in  the 
final  map  and  then,  after  plotting,  discarding  the  lines 
not  needed. 

In  summary,  it  appears  that  best  results  can  be 
obtained  with  the  program  by  specifying  a  larger  scale 
and  a  smaller  contour  interval  than  needed  for  the  final 
isohyetal  map.  After  plotting  the  isohyets  on  the 
computer  output  sheets,  only  those  contours  at  the 
prescribed  contour  interval  are  copied  and  reduced  in 
size  to  the  scale  needed. 


INPUT  AND  OUTPUT  OF  PROGRAM 


The  program  in  its  present  form  can  handle  an 
unlimited  number  of  events  for  which  isohyetal  maps  are 
required  for  a  given  watershed.  A  series  of  isohyetal 
maps  will  be  produced,  one  for  each  rainfall  event.  All 
maps  will  be  produced  at  the  same  scale,  and  each  will 
be  accompanied  by  tables  and  other  pertinent  informa- 
tion as  described  herein. 

The  structure  of  the  input  data  deck  and  the 
format  used  for  the  various  cards  are  described  below 
with  reference  to  figure  4.  The  deck  is  made  up  of  two 
main  parts:  Part  A  (fig.  4)  contains  information  about 
the  network  of  rain-gaging  stations;  Part  B  contains 
information  about  the  rainfall  events  to  be  plotted.  Part 
B  may  be  repeated  as  often  as  needed  to  include  all 
rainfall  events  for  which  isohyetal  maps  are  required. 
The  types  of  cards  used  for  the  input  data  deck  are 
described  below.  The  numbers  preceding  the  following 
paragraphs  correspond  to  the  numbers  of  the  various 
cards  shown  in  figure  4. 

1.  The  first  card  contains  information  about  the  water- 
shed and  the  rainfall  events  to  be  plotted.  The 
message  TITLE,  which  may  contain  up  to  80 
alphanumeric  characters,  will  be  printed  with  the 
isohyetal  map  and  the  tables  accompanying  it.  The 
card  format  is  (8A10). 


2.  The  second  card  controls  the  size  and  scale  of  the 
map  to  be  produced.  It  contains  the  following 
information: 

a.  YMIN,  YMAX,  XMIN,  XMAX-minimum  and 
maximum  values  of  Y  and  X,  respectively.  These 
values  determine  the  borders  of  the  map  pro- 
duced. 


Figure  4.-Structure  of  deck  of  cards  containing  input  data. 


b.  SCALE—  the  number  of  miles  represented  by  1.0 
inch  on  the  map.  This  value  determines  the  size  of 
the  map. 

c.  ROT-  a  variable  indicating  if  rotation  of  the  axes 
is  needed.  A  value  of  zero  or  a  blank  entry  will 
cause  no  rotation.  Any  value  other  than  zero  will 
cause  a  clockwise  rotation  of  the  axes  by  90°. 

d.  PRNT—  a  variable  controlling  the  printing  of  a 
table  of  interpolated  data  points.  A  zero  value  or  a 
blank  entry  will  suppress  the  printing  of  the  data. 
Any  other  value  will  cause  the  printing  of  a  table 
giving  the  coordinates  and  rainfall  depth  for  all 
interpolated  points. 

The  format  of  card  No.  2  is  (8F10.0). 

3.  The  deck  of  cards  following  card  No.  2  describes  the 
rain-gaging  station  locations  and  the  definitions  of 
the  lines  along  which  interpolations  will  be  made. 
One  card  containing  the  following  information  is 
provided  for  each  station. 

a.  NRG— the  identifying  number  of  the  rain-gaging 
station. 

b.  AX,  AY— the  X  and  Y  coordinates  of  the  rain- 
gaging  stations. 

c.  NA,  NB,  NC,  ND— the  identifying  numbers  of  the 
adjacent  four  stations  defining  lines  in  the  interpo- 
lation network.  If  fewer  than  four  station  numbers 
are  given,  interpolation  will  be  made  only  to 
numbered  stations  as  provided. 

The  format  of  these  cards  is  (110,  2F10.0,  415). 

4.  Tlie  card  following  the  station-description  cards  is  a 
blank  card  signifying  the  end  of  these  cards. 

5.  The  first  card  of  part  B  of  the  input  deck  is  a  card 
identifying  the  individual  events  for  which  isohyetal 
maps  will  be  prepared.  One  to  seven  events  may  be 
included,  with  each  event  being  identified  by  a 
10-character  alphanumeric  word,  EVENT.  The  event 
identifiers  will  be  printed  with  the  isohyetal  maps 
and  tables.  The  format  of  the  card  is  (7A10). 

6.  The  second  card  of  part  B  specifies  the  contour 
interval,  DELH,  to  be  used  for  computing  the 
isohyetal  points.  The  same  contour  interval  will  be 
used  for  all  events  listed  on  the  previous  card.  The 
format  of  the  card  is  (FIO.O). 


7.  The  deck  of  cards  following  the  contour  interval  card 
contains  rainfall  depth  data  for  the  one  to  seven 
events  for  all  rain-gaging  stations.  Each  card  has  the 
identifying  number  of  one  station,  NGA,  and  the 
corresponding  rainfall  depth,  RDP,  for  the  events  for 
which  maps  are  to  be  prepared.  Each  station  must 
have  a  record  for  each  of  the  events  since  a  blank 
entry  is  interpreted  as  zero  depth  of  rainfall.  Stations 
for  which  no  cards  are  provided  will  be  ignored.  The 
format  of  cards  in  this  deck  is  (110,  7F10.0). 

8.  The  last  card  is  a  blank  card  signifying  the  end  of  the 
depths-of-rainfall  cards. 

The  number  of  rainfall  events  included  in  part  B  of 
the  input  deck  may  be  less  than  the  seven  events  for 
which  space  is  provided.  The  events  included  must,  in 
this  case,  occupy  the  corresponding  number  of  positions 
on  card  Nos.  6  and  7,  starting  from  the  left-hand  side  of 
the  cards.  A  repetition  of  part  B  of  the  input  deck  may 
follow  the  original  deck  even  if  the  number  of  events  in 
the  original  deck  was  less  than  seven. 

The  AlO  format  of  card  Nos.  1  and  5,  used  to 
supply  identification  information,  is  the  largest  allowed 
on  the  computer  available.  In  adapting  the  program  to 
another  computer,  these  formats  may  have  to  be 
changed  to  agree  with  the  capacity  of  the  computer 
used. 

The  output  of  the  program  is  self-explanatory.  It 
consists  of  the  following  sections,  each  starting  at  a  new 
page  of  output. 

1.  The  first  section  contains  information  about  the 
identification  of  the  watershed,  boundaries  specified 
for  the  map,  and  the  scale  used. 

2.  The  identifications  of  the  watershed  and  of  the  event 
considered,  and  a  list  of  comments,  if  any,  that  were 
produced  during  the  first  division  of  the  program  are 
printed  in  the  second  section. 

3.  The  third  section  lists  a  table  of  the  original  data  used 
for  the  interpolation,  giving  coordinates  and  depth  of 
rainfall  for  each  station.  Also  included  in  this  section 
is  a  statement  about  the  contour  interval  and  the 
number  of  rain-gaging  stations  used. 

4.  Section  4  contains  a  table  of  interpolated  data,  giving 
coordinates  and  depth  of  rainfall  for  all  isohyetal 
points  interpolated  by  the  program.  The  points  are 
arranged  by  decreasing  value  of  the  Y  coordinate  or 
increasing  value  of  the  X  coordinate.  The  table  is 
followed  by  a  statement  of  the  number  of  interpo- 
lated points.  The  table  of  interpolated  data  will  be 


suppressed  unless  the  print  option  is  chosen.  If  the 
table  is  not  printed,  the  statement  about  the  number 
of  interpolated  points  will  be  printed  with  section  3. 

Section  5  is  a  map  of  isohyetal  points  represented  by 
letters  of  the  alphabet  and  stations  represented  by 
asterisks,  all  plotted  at  their  locations  according  to 
the  scale  selected.  The  map  is  followed  by  identifica- 


tion statements  and  notes  concerning  the  scale  of  the 
map  and  the  isohyetal  interval. 

6.  An  "end  of  data"  statement  is  printed  in  section  6  if 
all  data  have  been  processed  and  the  exit  is  normal.  If 
the  program  is  stopped  because  of  duplicate  data,  a 
statement  to  this  effect  will  be  printed. 


APPLICATIONS  OF  THE  PROGRAM 


The  use  of  the  computer  program  will  be  illus- 
trated with  data  for  the  storm  of  July  28,  1966,  on  the 
Walnut  Gulch  watershed.  The  area  for  which  the 
isohyetal  map  will  be  prepared  is  part  of  the  main 
Walnut  Gulch  watershed,  contributing  to  measuring 
flume  No.  11  (watershed  No.  63.011).  The  boundaries 
of  the  subwatershed  and  the  locations  of  rain-gaging 
stations  on  and  near  the  boundaries  of  the  watershed  are 
shown  in  figure  5. 

The  first  step  in  using  the  program  involves  joining 
adjacent  rain-gaging  stations  to  produce  the  interpola- 
tion network.  The  network  for  the  stations  included  in 
figure  5  is  shown  in  figure  6.  The  first  part  of  the  input 
data  deck  (part  A  of  fig.  3)  is  prepared  from  information 


in  figures  5  and  6.  It  includes  map  boundaries,  rain- 
gaging  stations,  coordinates  from  figure  5,  and  defini- 
tions of  interpolation  lines  from  figure  6.  Part  B  of  the 
input  data  deck,  containing  rainfall  data  for  the  evening 
storm  of  July  28,  1966,  is  prepared  from  the  precipita- 
tion records  for  the  rain-gaging  stations.  The  two  parts 
of  the  input  data  deck  are  Usted  in  table  1 . 

The  output  of  the  computer  program  is  shown  in 
table  2  and  in  figure  7.  Table  2  gives  the  information 
printed  as  secfions  1  and  3  of  the  output,  and  figure  7  is 
a  reproduction  of  the  map  of  isohyetal  points  produced 
in  section  5  of  the  computer  output.  Figure  8  shows  the 
same  map  with  the  isohyetal  lines  drawn  through  the 
isohyetal  points.  The  final  isohyetal  map  for  watershed 
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Figure  5. -Coordinate  system,  rain-gaging  stations  and  boundaries  for  watershed  63.01 1. 
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TABLE  1. -Listing  of  input  data  deck 


Input 
deck  part 

Card 
type 

Listing  of  Cards 

A 

1 

Walnut  Gulch  Watershed  Above  Flume  No. 

11 

2 

10.0             13.0               18.6                             23.3 

0.40 

3 

38             18.965            11.820            43 

43             19.200            12.335            51              90 

89 

44             19.195             10.540            38 

50             20.180            12.645            43              89 

88 

55 

51             19.565            11.360            38              44 

57 

52 

52             20.210            10.715            90             44 

57 

72 

54             21.290            12.020            50              88 

56 

55             22.315            11.915            54              91 

60 

56             21.225            10.980            88              52 

72 

57             20.080            10.185            72              44 

60             22.140            11.035            56              72 

61            22.000            10.485            72             56 

60 

64             23.170            12.275            55              60 

50 

54 

65             23.190            11.150            60              64 

55 

72             21.015             10.190            90 

82            22.400            10.540           60             61 

65 

88            20.780            11.705            90             89 

52 

89             20.100            12.085            38              51 

54 

90             20.170            11.465            51              89 

56 

4 

91            21.840            11.540           60             54 

88 

56 

B 

5 

July  28  66 

6 

0.1 

7 

38             1.28 

43             0.88 

44              1.34 

50              1.12 

51              1.08 

52             1.43 

54             1.88 

55             2.23 

56             1.68 

57             1.55 

60             1.93 

61              2.01 

64             2.15 

65             2.07 

72             2.12 

82             2.18 

88             1.54 

89             1.01 

90             1.19 

91              1.94 

8 

63.01 1,  prepared  on  the  basis  of  the  points  produced  by 
the  computer,  is  shown  in  figure  9.  The  area  represented 
in  figure  9  is  only  part  of  the  area  shown  in  figure  8. 

The  computer  time  needed  to  run  this  program  on 
the  CDC^  6400  computer,  using  Fortran  IV  language,  at 
the  University  of  Arizona  was  5.5  seconds  on  the  central 
processing  unit  and  20.4  seconds  on  the  peripheral 
processors.  These  times  include  the  compiling  of  the 
program,  which  took  3.9  seconds  on  the  central  process- 
ing unit  and  14.4  seconds  on  the  peripheral  equipment. 


Figure  6. -Interpolation  networK  lor  rain-gaging  stations  on  and 
near  watershed  63.01 1. 


Use  of  a  trade  name  does  not  imply  endorsement  by  the 
U.S.  Department  of  Agriculture.  It  is  used  only  for  the  purpose 
of  providing  specific  information  to  readers. 


TABLE  2. -Output  of  computer  program 


Walnut  Gulch  Watershed  Above  Flume  No.  1 1 

YMIN 

YMAX 

XMIN 

XMAS 

SCALE 

10.0000 

13.0000 

18.6000 

23.3000 

.4000 

Walnut  Gulch  Watershed  Above  Flume  No.  1 1 

July  28  66 

ORIGINAL  DATA  USED  FOR  INTERPOLATION 

RG 

X 

Y 

RAIN 

38 

18.9650 

11.8200 

1.2800 

43 

19.2000 

12.3350 

.8800 

44 

19.1950 

10.5400 

1.3400 

50 

20.1800 

12.6450 

1.1200 

51 

19.5650 

11.3600 

1.0800 

52 

20.2100 

10.7150 

1.4300 

54 

21.2900 

12.0200 

1.8800 

55 

22.3150 

11.9150 

2.2300 

56 

21.2250 

10.9800 

1.6800 

57 

20.0800 

10.1850 

1.5500 

60 

22.1400 

11.0350 

1.9300 

61 

22.0000 

10.4850 

2.0100 

64 

23.1700 

12.2750 

2.1500 

65 

23.1900 

11.1500 

2.0700 

72 

21.0150 

10.1900 

2.1200 

82 

22.4000 

10.5400 

2.1800 

88 

20.7800 

11.7050 

1.5400 

89 

20.1000 

12.0850 

1.0100 

90 

20.1700 

11.4650 

1.1900 

91 

21.8400 

11.5400 

1.9400 

INTERVAL  BETWEEN  ISOHYETS  =  .100  INCH 

THE  NUMBER  OF  RAIN-GAGING  STATIONS  USED  is 

20 

THE  NUMBER  OF  INTERPOLATED  POINTS  IS  189 
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SCALE    IN    MILES 


Figure  9.-Isohyetal  map  for  storm  of  July  28,  1966,  (1900)  on  watershed  63.01 1. 


The  net  time  for  the  preparation  of  the  map  is  thus  1 .6 
seconds  on  the  central  processing  unit  and  6.0  seconds 
on  the  peripheral  processors. 

The  following  information  is  given  as  another 
example  of  computer  time  needed  to  run  the  program. 
The  preparation  of  a  map  of  isohyetal  points  for  storms 
on  the  Walnut  Gulch  watershed  (63.001)  based  on  data 
from  90  rain-gaging  stations  took,  exclusive  of  compiling 
time,  about  5.6  seconds  per  storm  on  the  central 
processing  unit  and  about  9.6  seconds  per  storm  on  the 
peripheral  processors.  The  area  of  the  watershed  is  57.7 
square  miles.  The  scale  of  the  map  produced  was  0.75 
mile  per  inch,  and  the  contour  interval  was  0.2  inch.  The 


number    of  interpolated   points   was   about   400   per 
storm. 

The  computer  program  has  been  applied  to  a 
number  of  storms  on  a  few  watersheds  with  satisfactory 
results.  In  another  application  it  was  used  to  determine 
the  distribution  of  vegetation  and  soil  composition  on  a 
watershed.  The  percentage  of  soil  constituents  and  the 
density  of  given  species  of  vegetation  were  used  as  the 
variables  to  be  plotted,  rather  than  the  rainfall  depth. 
The  program  can  be  used  also  for  plotting  of  ground- 
water-table  contours,  lines  of  equal  depth  to  ground- 
water, equal-temperature  lines,  or  in  fact  maps  of  any 
variable  that  is  a  function  of  location. 


CALCULATION  OF  VOLUME 


The  computer  program  for  plotting  isohyetal 
points  does  not  and  cannot  compute  the  average  depth 
of  rainfall  or  the  volume  of  rainfall  deposited  on  the  area 
represented  by  the  map.  The  reason  is  that  the  map 
produced  contains  only  isolated  isohyetal  points, 
whereas  the  complete  isohyets  are  needed  for  the 
computation  of  the  average  depth  of  rainfall. 

It  is  generally  agreed  that  the  average  depth  of 
rainfall  derived  from  an  isohyetal  map,  prepared  by  the 
conventional  interpolation  methods,  is  very  nearly  equal 
to  the  average  depth  of  rainfall  computed  by  Thiessen 
weights.  If  the  isohyetal  maps  are  to  be  used  for 
computing  the  average  depth  of  rainfall,  it  is  necessary 


to  read  new  input  data  from  the  maps  for  each  storm. 
For  computation  based  on  Thiessen  weights,  all  that  is 
needed  is  a  simple  program  that  multiplies  the  observed 
depth  of  rainfall  by  Thiessen  weights  computed  once  for 
the  network  of  rain-gaging  stations. 

A  computer  program  for  computation  of  Thiessen 
weights  has  been  described  elsewhere.^  Copies  of  the 
deck  of  cards  of  this  program  are  also  available  from  the 
Southwest  Watershed  Research  Center,  Southwest 
Branch,  Agricultural  Research  Service,  Tucson,  Ariz. 
85705. 


Diskin,  M.  H.  "On  the  computer  evaluation  of  Thiessen 
Weights,"  Jour,  of  Hydrology   11(1):   135-144,  1970. 
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CONCLUSION 


A  computer  program  has  been  prepared  for  com- 
puting the  locations  of  isohyetal  points  and  their 
plotting  as  part  of  the  printed  output  of  a  digital 
computer.  The  locations  of  the  points  are  found  by 
linear  interpolation  along  lines  of  an  interpolation 
network  defined  by  the  input  data.  The  points  are 
plotted  as  alphabet  letters  printed  at  the  computed 
locations  for  the  map  scale  specified.  Each  letter  signifies 
a  different  isohyet  according  to  the  contour  interval 
selected. 

The  map  borders  have  X  and  Y  coordinate  scales 
marked  at  regular  intervals.  The  map  is  accompanied  by 
explanatory  notes  and  tables  of  data.  An  option  is 
available  in  the  program  for  rotating  the  axes  by  90° 


before  printing,  if  this  is  advantageous  for  the  shape  of 
the  watershed  considered.  The  locations  of  the  rain- 
gaging  stations  are  also  included  in  the  map. 

The  program  reduces  the  manual  work  involved  in 
the  preparation  of  isohyetal  maps.  It  produces  results 
comparable  to  or  better  than  the  manually  prepared 
maps.  The  computer  time  needed  for  running  the 
program  is  relatively  short.  The  number  of  maps 
prepared  in  one  setting  is  practically  unlimited.  In 
addition  to  the  preparation  of  isohyetal  maps,  the 
program  may  be  used  for  the  preparation  of  maps  of 
ground-water-table  contours  and  other  contour  maps 
useful  for  hydrologjcal  investigations. 
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APPENDIX 

PROGRAM  CoNPLT  ( InpUT, OUTPUT* TAPf 1 =INPU r ) 

c 

C  A  PROGRAM  TO  INTERPOLATE  ISOHveiTAL  POINTS  BETWEEN  RAiNbA(itS 

Ci         AND  To  PLOT  THE  INTERPOLATED  POINTS  aT  a  GIVEN  SCALE. 

C 

C  PROGRAM  CONPLT  WAS  WRITTEN  BY  'Uh.UlSKiN, 

C 

C  THE  PROGRAM  IS  HASEU  ON  PROGRAMS  CONTOUk  ANO  MAP  wkITIEN  oY 

Cl         P.R.OAVIS.  UNIVERSITY  OF  ARIZONA. 

C 

REAL  AX (200)  f AY  (2c 0)  »HR(20U) fX (3000)  »Y(J0O0)  »SCL(12)  » 
iHOuOO)  ♦XS(  3000)  •VQOnOO)  ,HS(  3000)  tRDP(?00*l0) 
INTEGER  N|jG(2no)  ♦ma  (  2no  )  .NB  (200)  tNC  (200  )  ,ND(200)  9J(U^) 
INTEGER  NGA(200)  tL-ST(loOO) 
DIMENSION  TITLE(8),  EVENT(8).  ALPHA ( 121 ) •rtET A ( 121 ) 

READ  115«  TITLE 

READ  1C1»  YMlN,Yl^AX»XMlN,xHAX,SCA|.E»ROT»PRNr 

PRINT  100 

PRINT  115.  TITLE 

PRINT  117 

PRINT  114 

PRINT  116t  YMIN,YMAXtXMlN,XMAX»SCALE 

PRINT  117 

DELXasCALE/10,0 

0£LY=SCALE/6.n 

DO  66  1=1.1000 

LST(I)=0 
66    CONTINUE 

00  37  1=1,200 

HR(I)so.O 

NRG(I)=0 

NGA(I)=0 
37    CONTINUE 

PTX=0.0 

DO  51  1=1.200 

READ  102»  NRG(I)  ./^-Xd)  .AYd)  »NA(I)  .NBd)  ,NC(I)  tNU(I) 

IF{NRG(I) .EQ.o)  GO  TO  59 

MRsI 

JaNRG(I) 

IF(LST(J)  ,NE.f>)  60  TO  80 

LST(J)=I 
51  CONTINUE 
59    IF(ROT.EQ,o.o)  GO  TO  60 

DO  75  1  =  1. -vIR 

TEMP=AX(I) 

AX(I)»AY(I) 

AY(I)=-TEMP 
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75 

CONTINUE 

TemPsXMIN 

TEMA'XMAX 

XMINsYMIN 

XMAX=YMAX 

ymins-tema 

YMAXa-TEMP 

6^ 

READ  U5t  (EVENT(J)  »J  =  1  »7) 

IF(EOF»n  73»74 

74 

READ  lOlf  DELH 

NAV  =  2999-vlR 

00  34  I=l«200 

LST(I)=IABS(LST(I)) 

34 

CONTINUE 

DO  68  1=1.200 

HEAD  105»NGA(I) f (PnP(l»J) .J=lf7) 

IF(NGA(I) .EQ,0)GO  TO  35 

JRsNGAd) 

IF(LST(JR) .LT.O)  GO  TO  81 

LST(JR)=-LST(JR) 

MGAsI 

68 

CONTINUE 

35 

DO  50  NEV=l,7 

IF  (NEV,E0,4)  GO  TO  50 

PRINT  100 

PRINT  141»TITI.E.EVENT(NEV) 

PRINT  117 

Do  70  1=1. MGA 

NRaNGA(I) 

RFsROPd.NEV) 

IF  (I.EQ.l)  HMX  =  RF 

lF{RF,GT,HMx)  HMX  =  RF 

IF(LST(NR) .EQ.O)  PRINT  118.NR 

IF(LST(NR)  ,EU.O)  f^-0  TO  70 

JRa-LST(NR) 

HR(JR)=RF 

70 

CONTINUE 

PRINT  117 

IF(HMX,EQ,0»0)  GO  TO  60 

HST  =  HMX  ♦  l00«0*nELH 

ALM«DELH«0,001 

NPb^ 

DO  54  1=1. MR 

JB«NRG(I) 

IF(LST(JB),GE,0)  PRINT  12n.  J8 

IF(LST(JB) .GE.O)  GO  To  54 

J6(l)s  NA(I) 

JG{2)x  NB(I) 

JG(3)a  NC(I) 

J6(4>a  ND(I) 

DO  55  JBI.4 
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IF  (JG(J) ,EQ,0)  00  TO  S5 

JMsJG( J) 

IF(LST(JM) ,GE.O)  PHINT  12l»JBfJM 

IF(LST(JM)  .GE.O)  ^-0  Tn  55 

JAs-LST(JM) 

DIFX=AX(JA)-AX(I) 

DIFHsHR(Ja)-HR(I) 

DIFYSAY(JA)-AY(I) 

DABHaABS(DlFH) 

IF(nABH,LT.ALM)  GO  TO  55 

IF(HR(I) .GrtHRC JA) )  GO  TO  61 

XCL=HR(I)/OELH-0.? 

XCVsHR(JA)/DELH*0,? 

Go  TO  62 
61    XCLeHR(JA)/DELH-0,? 

XCV=HR( I)/DELH*0.? 
6?    NCL=XCU*1.0 

NCVsXCV 

IF(NCL.GT,NCV)  GO  TO  55 

DO  56  K=NCL»NCV 

NPsNP*! 

H(NP)sDELh«FLOAT(K) 

X(NP)sAX(I)*(H(NP)-HR(I) )»DlFX/UlFH 

Y(NP)=AY(I)*(H(NP)-HR(L) )»DlFY/DlFH 

IF  (AX(I),EQtAX(JA))  GO  TO  79 

IF  (X(NP)  .LE.AX(Ja)  .ANO,X(NP)  .GE.AXd)  )  GO  TO  49 

IF  (X(NP)  •GE.AX(JA)  .ANO.X(NP)  •LE.AXd)  )  GO  TO  49 

XID=ABSU(NP)-AX(T)  ) 

XJD=ABS(X(NP)-AX(JA) > 

XOL=2tO*DfcLX 

XDMaABS(0,C3*DlFX) 

GO  TO  67 
79    IF(Y(NP)  .LE.AYlJA)  .AND.Y(NIP)  tGE.AYd)  )  GO  TO  49 

IF(Y{NP) .GE.AY(JA) .AND»Y(NP) tLE.AY(I) )  GO  TO  49 

XIDsABS(Y(NP)-aY(I)) 

XJD=ABS(Y(NP)-aY(Ja)) 

XDL=2.0»DELY 

XDM«ABS(0,03*OIFY) 
67    IF(xJD,LT.XlO)XlD=xJn 

IF(XDL.LT,XOM)  XOL=XOM 

IF(XI0,GT,XDL)NP=NP-I 

49    CONTINUE 

IF(NP,EQ«NAV)  GO  TO  64 

56  CONTINUE 
55  CONTINUE 
54  CONTINUE 
45    PRINT  100 

PRINT  141»TITLE»EVENT(NEV) 

PRINT  117 

PRINT  128 

PRINT  126 

Mao 

00  46  1=1. MR 
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JaNRG(I) 

IF  (LST( J) .GE.O)  GO  To  46 

NP  =  NP  ♦  1 

X(NP)  =  AX(I) 

Y(NP)  s  AY(I) 

H(NP)  =  HST 

TEMP=-flY(l) 

IF  (ROT.E'J»0.0)  PRINT  110,  NRG(I)  f  AX(I)  f  AY(I)  ♦HR(I) 
IF  (ROT.Ne.O.O)  PWiriT  110,  NRG(I),  TEMP, /\X  (  I )  ,HR  (I ) 
46    CONTINUE 
NPsNP^l 

X(N1P)=XMIN*0ELX«5,-' 
Y(NP)sYMlM-DELr»5,o 
H(NP)=H5T 
PRINT  117 
PRINT  123»  DELH 
PRINT  117 
PRINT  135»M 

IF  (PRNT.EQ.O)  GO  To  33 
PRINT  100 

PRINT  141,TlTLt,tVF  -ITCNEV) 
PRINT  117 
PRINT  129 
PRINT  127 

33    M»o 

DO  57  J=1»NP 

YS(J)=YMIN-DELY»8.n 

DO  58  Isl,np 

IF(Y(I) .LT.YS(J) )  GO  TO  5R 

YS(J)=Y(I) 

NSsI 
58    CONTINUE 

Y(NS)sYMlN-DELY»l,'%'J 

XS(J)=X(NS) 

HS(J)=H(NS) 

IF  (HS(J) .EQ.HST)  GO  TO  57 

MsM*l 

IF  (PRNT.EQtO.O)  GO  To  57 

TfMPs-YS(J) 

IF  (ROT.Eino)  PRINT  110»M,XS( J) fYsC J) tHSCJ) 

IF  (ROT.NE.O)  PRirJT  1 1  0  ♦M,TeMP,XS  (  J)  fhS  ( J) 
57    CONTINUE 

PRINT  117 

PRINT  132«M 

00  31  I=l»121 

BETA(I)slH- 

31  CONTINUE 

DO  32  Isl»121»lO 
BETA(I)slH* 

32  CONTINUE 
WESTsXMiN 


28 


30 


38 


36 
42 

40 


?9 
41 

63 


EAST»WEST*120.0»UELX 

SCL(1)=WEST 

DO  30  I=2»12 

SCL(I)=SCL(I-l)*l;^.0«nELX 

CONTINUE 

PRINT  151, SCL 

PRLsYMAX-  •5«DELY 

pnl=ymax 

IF(ROT,NE»0.0)  PNL=-PNL 

PRINT  137»  PNL.BETA 

LbC 

PRL=PRL-DELY 

PNL5PRL*0,5»OELY 

IF  (ROT. NE. 0.0)  PNL=-PNL 

DO  36  I=lfl21 

ALPHA(I)=1H 

CONTINUE 

IF(YMIN.GE.PRL)  ^50  TO  48 

L=L*1 

IF(YS(L)  .LT.PRL)  (50  TO  43 

IF(XS(L) .LE. EAST. AND. XS(L) .GE. WEST)  GO  TO  40 

60  TO  42 

X I« ( XS (L) -WEST )/OrLX* 0,5 

I=XI*1 

IF(ALPHA(I).NE.1H   )  I=I*i 

IF(ALPHA(I) ,NE.1H   )  1=1-2 

IF(I.LT.l)  1=1 

IF(I.GT.121)  1=121 

IF  (HS(L) .EQ.  HST)  GO  TO  63 

XDIP^=HS(L)/0ELH 

KDIFaXniF 

RmD=XDIF-M)IF 

IF  (RMn.GT.0.5)KniF  =  Kr)IF*l 

IF  (KDIF,EO.O)GO  TO  26 

IF(KDIF.LE.26)  GO  TO  41 

IMeM00{KDIF.26) 

IF(IM,NE.C)  GO  TO  -^9 

GO  TO  26 

NK=K0IF/26 

K0IF=K0IF-NK»26 

GO  TO  (l»2»3»4»5.6.7»rt»9,lO,11.12,  I3»14,i^,l6,l7tl'^» 

il9,20»21»22»2i 

lt24,25,26) ^KDIF 

ALPHA{I)=1H» 

GO  TO  42 

ALPHA(I)=lHA 

GO  TO  42 

ALPHA(I)=lHR 

60  TO  42 

ALPHA{I)=IHC 

GO  TO  42 

ALPHA(I)=1HD 

GO  TO  42 
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5  ALPHA{I)=1HE 
GO  TO  k?. 

6  ALPHA(I)=1HK 
GO  TO  4? 

7  ALPHA(I)=lHG 
GO  TO  43 

8  ALPHA(I)=lHH 
GO  TO  42 

9  ALPHA(I)=1HI 
GO  TO  4? 

10  ALPHA(I)=1HJ 
GO  TO  42 

11  ALPHA(I)=]HK 
GO  TO  4? 

12  ALPHA(I)=1HL 
GO  TO  t^Z 

13  ALPHA(I)=1HM 
GO  TO  42 

14  ALPHA(I)=1HN 
GO  TO  42 

15  ALPHA(I)=1H0 

60  TO  42 

16  ALPHA(I)slHP 

GO  TO  42 

17  ALPHA(I)=1HQ 
GO  TO  42 

18  ALPHA(I)=1HR 

GO  TO  42 

19  AUPHAa)=lHS 
GO  TO  42 

20  ALPHA(I)31HT 
GO  TO  42 

21  ALPHA(I)=1HU 
GO  TO  42 

22  ALPHA(I)=1HV 
GO  TO  42 

23  ALPHA(I)=1HW 
GO  TO  42 

24  ALPHA(I)=1HX 
GO  TO  42 

25  ALPHA(I)siHY 
GO  TO  42 

26  ALPHA(I)slHZ 
GO  TO  42 

43    PRINT  145. PNLf ALPHA 

L=L-1 
GO  TO  38 
48    PRINT  137.  PNL.  BETA 
PRINT  113.  SCL 
IF  (ROT.EQ.O.O)  PRINT  133 
IF  (ROT.NE.O.O)  PPlNiT  134 
PRINT  141,  TITLE. FVENT(NEV) 
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IF  (PTX.GT.O.C)  PkINT  12^5 

PRINT  in 

PRINT  122.  SCALE»1'F:LH 

PRINT  117 

ASTP=1H» 

PPTNT  130.ASTR 

PHI^T  131 

iF(EAST.Gt.XMAX)  GO  TO  50 

WESTsEAST 

GO  TO  p.a 
5?.    CONTINUE 

GO  TO  6  0 
73    PRINT  142 

GO  TO  27 
64    PRINT  125 

PTXsUo 
GO  TO  45 

80  PRINT  124. J 
GO  TO  82 

81  PRINT  119.  NR 

82  PRINT  140 
27    STOP 

100  F0RMAT(»1») 

101  FORMATCeFlO.O) 

102  FORMAT(IlO.2F10.n,4l5) 
105   FORMAT (110 f7F 10.0) 

110   FURMAT(I12»3F15.4) 

113  F0PMAT(/»7X,12Fl(i,3.) 

114  FORMAT (7X,4HyMlN»8X.4HYMAX.8X,4HXM IN, rtX.^HXMAX. /X.SHSCALfc.) 

115  F0RMAT(8Alr^) 

116  F0RMAT(6F12.4) 

117  F0RMAT(//) 

118  FORMAT  (5X.21HRAINFALL  FOR  OAGF.  NO.»I^» 

1  iflH   NOT  USED.  COORDINATES  OF  GAGF.  NOT  GIVEN) 

119  Format (5x,5oHMoRE  than  one  rainfall  value  is  given  for  gage  no., 

115) 

120  F0RMAT(5X,35HN0  RaINFaLL  QATA  GIV^N  FOR  GAGE  NO. ,15) 

121  F0RMAT(5X,?7HINTERP0LaTI0N  FROM  GAGE  N0.,I5,13H   TO  GaGE  nO..I5. 
1  23H   CANNOT  BE  CARRIED  OJT) 

122  F0RMAT(5X,14HSCALF  OF  MAP   ,Fl0.4,l6H   MILES  =  1  INCH,30X, 
1  27HINTERVAL  RET^fFN  ISOHYETS  =,Fq,3»6H   INCH) 

123  F0RMAT(5X,h7HINTERVAL  BETWEEN  ISOhYETS  =,F9.3,6H   INCH) 

124  FORMAT (5X,23HC00RDINATES  OF  GAGE  MO, 15, 
1  26H   ARE  GIVEN  MORE  THAN  ONCE) 

125  FORMAT(l0X,46HNUMnErK  OF  INTERPOLATED  POINTS  EXCEEDS  CAPACITY) 

126  FORMAT ( 1  OX, 2HRG,l?X.lHX,l4X,lHY,l0X,4HRA IN/) 

127  FORMAT (10X,2HPT,l?X.lHXfl4X,lHY,lnX.4HRA IN/) 

128  FORMAT(10X,36HORIGINAL  DATA  USED  FOR  INTFRPOLaTIOn//) 

129  FORMAT(IOX.12HINTFRI^OLATED  DATA//) 

130  F0rMAT(5X,32HrAIN  GaGE  LOCATION  DENOTED  QY    ,1A1  ) 

131  F0RMAT(5X,36HIS0HYFTAL  POINTS  DENOTED  BY  LETTERS  ) 

132  FORMAT (5X,38HTHE  NUMBER  OF  INTERPOLATED  POINTS  IS   ,I6) 
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133  Format  (/9ox»13hx  axis  miles/) 

134  FORMAT  (/9CX»13HY  AXIS   MfLES/) 

135  F0RMAT(5X,44HTHE  niUM8ER  Op  RAIN  G&GING  STATIONS  USEO  IS   .IS) 
137  F0RMAT(1X,F8,3,3X,1HI,121AI,1HI)  •  ^ '  i^  ^5  ^btU  IS   ,Ib) 

140  FORMAT (5X,32H  DUPLICATE  DaTa,  PROGRAM  SFoPPEU) 

141  FORMAT(5X,8A10,/,SX,1A10) 

142  FORMATdOXtllHENO  OF  DATA) 

145  F0RMAT(1X,F8.3,3X,HI,121A1,1HI) 

151  FORMAT (lHl,7X,12Fin. 3,/) 
END 
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